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Head-on infall of two compact objects with arbitrary mass ratio is investigated using the multipolar 
post-Minkowskian approximation method. At the third post-Newtonian order the energy flux, in 
addition to the instantaneous contributions, also includes hereditary contributions consisting of the 
gravitational-wave tails, tails-of-tails and the tail-squared terms. The results are given both for 
infall from infinity and also for infall from a finite distance. These analytical expressions should 
be useful for the comparison with the high accuracy numerical relativity results within the limit in 
which post-Newtonian approximations are valid. 
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The spiraling coalescence of two compact objects (black holes or neutron stars) moving about one another in an 
orbit, forms a prominent class of sources of gravitational radiation Such sources of gravitational waves (GW), 
especially in their late stages of evolution are prime targets for gravitational wave detectors such as LIGO Q and 
^ | Virgo [H]. The evolution of the binary systems composed of two compact objects involves three stages of evolution; the 
I 1 early inspiral, late inspiral and merger and the final ringdown. Detection of gravitational radiation from such systems 
vTrV ^ the gravitational wave detectors depends strongly on the theoretical inputs, which will involve computation of 
i ' the waveform of the signal for all the three phases to very high post-Newtonian (PN) order to detect and infer the 
characteristics of the sources of GWs, using matched filtering techniques Q. Even though head-on collision of two 
black holes has only a small astrophysical possibility, it provides the simplest possible situation to study the two-body 
problem of general relativity and has been studied since it provides an excellent theoretical platform for comparing the 
validity of various analytical and numerical approaches towards solving Einstein's equations in dynamical situations. 

One of the earliest attempts to solve the problem of head-on collision using a complete general relativistic approach 
was due to Davis ct al @. They discussed the emission of gravitational radiation due to the radial infall of a test 
particle in Schwarzschild spacetime from infinity, using Zerilli's equation for black- hole perturbations @. Because of 
the axial symmetry of the system, the problem simplifies considerably and yet retains the features of astrophysical 
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■ interest such as emission of gravitational radiation at infinity. In addition to this, head-on collision can be considered 
t—{ , as an approximation to the last stage of the inspiralling coalescence-when two objects merge together to form a 
!* ■ single object. The first attempt to solve the head-on collision of two equal mass black holes numerically was due 
. £^ \ to Smarr and Eppley 04^1 • This program has undergone substantial improvement in accuracy and reliability with 
advances in the understanding of numerical issues in the treatment of Einstein's equations and availability of better 
computing jToj . The head-on collision of two black holes with arbitrary mass ratio has been investigated numerically 
in [ll|, [l2[ and semianalytically (T^J. In a recent work [l4[ head-on collision of two equal mass, nonrotating black holes 
with ultrarclativistic speeds have been studied using numerical methods. The main result of this analysis is that in 
such a process (where the initial energy of the system is dominated by kinetic energy of black holes) the total amount 
of energy converted to gravitational waves is about 14% of the initial mass-energy for the system and corresponds to 
large luminosities of the order of 10~ 2 c 5 /G. Another study related to the collision of two equal mass, nonrotating 
black holes moving at ultrarelativistic speeds and with generic impact parameter [l5| suggests that such collisions can 
produce black holes rotating close to the Kerr limit and the energy radiated in such a process would be roughly 35% 
of the center-of-mass (CM) energy. 

Another approach which may be used to study the head-on collision of two compact objects is the PN approximation 
approach. Though PN methods are valid for arbitrary mass ratios, they eventually break down under situations like 
strong gravitational fields and high speeds. Simone, Poisson and Will (SPW) [3 investigated the problem of head- 
on infall and compared the PN approach with black-hole perturbation (BHP) theory They provided 2PN accurate 
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expression for the far-zone GW energy flux and showed, in particular, that the energy radiated during the infall is 
wellestimatcd by the quadrupole approximation combined with the exact test-body equations of motion (EOM) in 
Schwarzschild background. Also in a recent study [T3|, a hybrid method using both PN approximations and BHP 
theories has been used to study the head-on collision of two black holes and found that PN and BHP theories can 
explain the main features of gravitational radiation for head-on mergers. 

In this paper we investigate the problem of head-on infall using the multipolar post-Minkowskian (MPM) ap- 
proach [1 81 - 4231 ] and provide the complete 3PN accurate expression for the GW energy flux emitted during the radial 
infall of two compact objects towards each other. In addition to the simpler instantaneous part of the energy flux 
we also compute the more complex hereditary contributions up to 3PN order which involves the contributions due 
to tails, tails-of-tails and tail-squared terms. We discuss the head-on problem both for infall starting from rest at an 
initial finite separation (denoted case I) and similarly for infall starting from rest at infinite separation (denoted by 
case II). Instantaneous contributions at 2.5PN order and at 3PN order, computation of tails at 2.5PN, tail-of-tail and 
tail-squared terms at 3PN order are the new results of this paper. Our computations suggest that the total energy 
radiated in the process of head-on infall of two compact objects with equal masses is roughly about 0.0074% of the 
Arnowitt, Deser, and Misner (ADM) mass of the binary and the peak luminosities are typically less than of the order 
5 x 10~ 6 c 5 /G. Comparing our PN estimates with the numerical relativity results we can see that the PN estimates 
are smaller than the numerical results typically by a factor of 27 consistent with the expectation that a larger fraction 
of energy radiated indeed comes from the merger phase of the infall rather than from the early inspiral. 

This paper is organized in the following way. In Sec. |TT] we begin by providing the structure of the far-zone GW 
energy flux at 3PN order, relations connecting radiative multipole moments to source multipole moments and the 
decomposition of the expression for energy flux into instantaneous and hereditary contributions. Section [IIII lists the 
3PN EOM as well as the 3PN accurate expression for the center-of-mass energy in standard harmonic coordinates 
for the head-on case. In Sec. IIVI we give the expressions for the desired multipole moments at the PN order required 
for the computation of 3PN energy flux for head-on infall case. In Sec. [V] we first exhibit the instantaneous part 
of energy flux up to 3PN order in standard harmonic coordinates followed by the corresponding expressions in two 
alternative coordinates for possible comparison with numerical relativity results: modified harmonic (MH) and ADM. 
Section PVTl describes the computation of the hereditary part of the energy flux. Finally, in Sec. IVII1 we bring together 
the complete 3PN accurate expression for energy flux in ADM coordinates and the energy radiated during infall to 
some fixed radial coordinate. Section IVIIII contains a graphical display of the salient features and our conclusions. 
These results should be useful to compare and match to simulations using numerical methods in regimes where both 
treatments are expected to be the valid. The paper ends with a short appendix relating the expression for conserved 
energy in standard harmonic (SH) coordinates to that in ADM coordinates. 



II. THE FAR-ZONE GW ENERGY FLUX 



We start the discussion by writing the 3PN expression for far-zone GW energy flux in terms of the symmetric 
trace-free radiative multipole moments [U,[25|. The PN structure for GW energy flux reads as, 
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In the above expression Ul and Vl (where L = ii»2 • • • *z represents a multi-index composed of / spatial indices) 



are the mass-type and current-type radiative multipole moments respectively and U L and V^' l> denote their 
time derivatives. The moments appearing in the formula are functions of retarded time U = T — R/c in radiative 
coordinates. 

Equation (|2.1j) is the general formula for the computation of 3PN accurate energy flux for any general isolated 
source. In a recent paper [25l ] the complete third post-Newtonian energy flux has been computed for inspiralling 
compact binaries moving in quasi-elliptical orbits. In the present work we specialize to the case of head-on infall 
and compute the 3PN accurate far-zone GW energy flux emitted due to head-on infall of two compact objects with 
arbitrary mass ratio using the MPM approximation method. The radiative current-type moments (Vl) are related to 
the source current moments Jl whose expansion at each PN order contains the orbital angular momentum J which 
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vanishes in the head-on case. Thus the current-type moments Vl will not contribute to GW energy flux and for the 
head-on case, Eq. (|2.f [) reduces to the following form, 
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In the MPM formalism, the radiative moments U l and Vl are related to canonical moments Ml and Sl respectively 
and these canonical moments are in turn expressed in terms of source moments {II, Jl, Wl, Xl, Yl, Zl}- Since in the 
present work we only deal with head-on situation we would exclude terms involving current-type multipole moments 
from all our expressions for the reason stated above. It should be evident from the Eq. (|2.2[) that for the computation 
of 3PN accurate energy flux Uij is needed at 3PN order, is needed at 2PN order, E/y/.; with 1PN accuracy and 
Uijkim t° leading Newtonian accuracy. General expressions connecting Ul to source moments have been listed in p5| 
and we shall simply recall those expressions. For the 3PN accurate mass quadrupole we have 
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where the bracket <> surrounding indices denotes the symmetric trace- free projection. The Il& are the mass- type 
source moments (and 1^ denote their n th time derivatives), and W is the monopole corresponding to the gauge 
moment Wl which for our purpose needs to be known Newtonian accuracy. The quantity M appearing in the 
above expression is the ADM mass of the source. It should be evident from Eq. (|2.3p that radiative moments have 
two distinct contributions. The first referred to as the instantaneous contribution requires the knowledge of source 
multipole moments only at a given retarded time, U = T — R/c; where R is the distance of the source in radiative 
coordinates. The second one, referred to as the hereditary contribution, which is given by integrals over retarded 
time from to 00, depends on the dynamics of the system in its entire past history and requires the knowledge of 
source moments at all times before U. A closer look at the hereditary terms reveals two types of contributions, some 
with and some without the log factors. The integrals (with log factors) appearing at 1.5PN and 3PN order are called 
tail and tail-of-tail integrals respectively. The integral (without log factor) appearing at 2.5PN order is called the 
nonlinear memory integral. It is a time antiderivative and hence leads to an instantaneous term in the energy flux. 
The mass-type octupole moment Uijk which is needed at 2PN is related to the associated source moment as 



u ijk (u)=imu)+ 2 -^- r 



. ct\ 97 
1Og (27j + 60 



(2.4) 



For other radiative moments, UijM, and Uijhimi only the leading order accuracy in the relation between radiative and 
source moments is needed, so that 



U L {U) = I<i\lJ) 



0(3). 



(2.5) 



The constant r which provides a scale for the logarithmic term in the above expressions is an arbitrary constant. It 
enters the relation connecting retarded time U = T — R/c in radiative coordinates to retarded time u = t — r/c in 
harmonic coordinates (where r is the distance of the source in harmonic coordinates) . The relation between retarded 
time in radiative coordinates, and the one in harmonic coordinates reads as 
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Later in this paper we shall show that the presence of this constant r will not influence any physical result like 
far-zone GW energy flux. 
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We can now use the expressions for the radiative moments given by Eqs. (|2.3[) - (|2.5[) in Eq. (|2.2[) to obtain the 
3PN energy flux formula in terms of source moments. As discussed above, the presence of two distinct contributions 
(instantaneous and hereditary) leads to a natural decomposition of the 3PN energy flux into two pieces and the 
complete flux can be written as a sum of the two distinct types of contributions as 
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where the instantaneous contribution 1 to the energy flux is given by 
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The hereditary contribution comprises of three parts, 

-7"hercd = J~ bail + -F tail(tail) + -F (tail) 2 

where the quadratic-order (proportional to G 2 ) tails are given by 
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and the cubic-order tails (proportional to G 3 ) by 
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Here one should note that the general formulae for energy flux include some contributions from current-type moments 
as well (see (2f| H3|) but these vanish for the head-on case. Further, it should be noted that Eqs. (|2.3p - (|2.4[) and 
thus Eqs. (|2.10j) - (|2.11j) show an intermediate dependence on the arbitrary length scale tq which should eventually 
cancel from all physical quantities. Such a cancellation of the scale tq from all physical quantities occurs naturally 
in the MPM formalism and has been explicitly shown for sources such as binary systems moving in circular [20l ] and 
elliptical orbits J2(| [27||- This is facilitated because an explicit computation of the hereditary integrals is possible 
since the integral over the complete past in the adiabatic approximation reduces to an integral over the current 
(noninspiralling) orbit which can then be computed making explicit use of the periodicity features in the motion. For 
a head-on situation, on the other hand, the absence of periodicity prevents the straightforward extension of the above 
method. A more first-principle treatment is called for based on the observation that since most of the ro- dependence 
comes from our definition [Eq. (|2 .6[) ] of a radiative coordinate system, it can be tracked and isolated by inserting U 
as given by Eq. (|2.6p back in Eq. (|2.3p - (|2.4[) . Upon doing so we get expressions for radiative moments in harmonic 
coordinates (t, r), which read 
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In the 1.5PN term, the ro dependence is more trivial and disappears with the change from radiative to harmonic 
coordinates. At 3PN order, however, there still remains a nontrivial r dependent term. However, the quadrupolc 
mass moment ly also depends on the constant ro at 3PN order [see Eq. (|4.2[) ] and when one takes those dependences 
into account we will check that the 3PN radiative moment Uij is indeed independent of ro. Using the expressions 
for radiative moments given by Eqs. (|2.12[) and (|2.13j) we can rewrite explicitly the different hereditary contributions 
given by Eqs. (|2.9j) . (|2.10|) . and (|2.1ip . One finds the quadratic-ordcr(proportional to G 2 ) tails are given by 
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We shall come back to the discussion of the hereditary terms in detail in Sec. I VII where we shall compute their 
contributions to far-zone energy flux. 



III. THE EQUATIONS OF MOTION AND THE CONSERVED ENERGY FOR HEAD-ON COLLISION 

A. Standard Harmonic Coordinate System 

By standard harmonic coordinates we refer to the coordinate system that has been used in previous works [l8l [28j . 
Since the head-on collision problem has only one direction of motion one can write relevant equations for the head-on 
case by imposing the restrictions 

x = zn, v = in, r = z, v = r = z, (3-1) 

on the corresponding expression for general orbits in terms of z and i, where z is the separation between the two 
objects at a given time and i is the first time derivative of z, giving (coordinate) speed with which they are moving 
with respect to each other at that instant. 

The computation of the energy flux involves time derivatives of source multipole moments which in turn will require 
the knowledge of equations of motion at appropriate PN order. Computing the 3PN accurate energy flux requires the 
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3PN accurate equations of motion |28| . 3PN accurate equations of motion in the CM frame associated with standard 
harmonic coordinate system, for compact objects moving in generic orbits, are given in [25l. [28j. 

Since we will discuss the results in other coordinate systems like the modified harmonic coordinates and the ADM 
coordinates in the subse quen t sections, we will provide for a more general discussion of the 2.5PN terms along the 
lines of [2!| based on [30ll3l| (see also [IH). In addition to the contact transformation involving "conservative" orders 
up to 3PN required to go from the SH to the modified harmonic and ADM coordinates (involving even order 2PN 
and 3PN terms) there still remains the possible change of gauge in the radiation reaction (dissipativc) terms at order 
2.5PN. Recall, that in the SH coordinate system the lowest-order dissipativc part of the equations of motion, i.e. the 
2.5PN acceleration term, is given by (with boldface letters indicating ordinary three-dimensional vectors) 
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One may however prefer to employ alternative radiation gauges and a convenient characterization at 2.5PN order has 
been investigated earlier in [30, |3l[ ■ Following this work, the most general form of the 2.5PN term in the relative 
acceleration is specified by the two-parameter family written as, 
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The general 2.5PN gauge is parametrized by the two numerical constants a and /3. The SH (and modified harmonic) 
gauge in which the acceleration is given by (|3.2[) corresponds to the choice a = — 1 and (3 = 0; the ADM gauge 
corresponds to a = 5/3 and (3 = 3, in which case the 2.5PN acceleration becomes [33| 



ADM 
a 2.5PN 



G 2 m 2 u 



12 v 2 



Gm 



15r z 



11 



1 Gm 
3~ 



Sr 2 



(3.4) 



By imposing the restrictions given by Eq. (|3.1[) we can write the equations of motion (or acceleration) in terms of 
the variables z and i for head-on situation as a 1 = an 1 , where, 
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The general expression for center-of-mass energy E associated with standard harmonic coordinate system is given 
in [28[ and the corresponding expression for head-on situation can be obtained by imposing restrictions given by 
Eq. (jXTj) . Thus we have 
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To study the head-on infall we consider two different situations, following [lfj. In the first case (we will call it case I) 
we assume that the radial infall proceeds from rest at a finite initial separation whereas in the second case (case II) , 
we assume the objects start falling towards each other from rest at infinite separation. 

1. Case I: Infall from a finite distance 

Let us suppose the two objects initially separated by the distance z, start falling radially towards each other from 
the rest i.e. z {z{) = 0. Hence the center-of-mass energy E in standard harmonic coordinates at z% will be 
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Gm/ziC . Equating, Eqs. (|3.6|) and (|3.7|> . the resultant expression can be inverted for z. 
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where 7 = Gm/zcf is the PN parameter and s = z/zi < 1. 



2. Case II: Infall from infinity 



We can view case II as a limiting case of case I and the expression for z can be obtained by inserting s = z/zi back 
in Eq. (|3.8[) and taking the limit when Zj — ¥ 00. For 2 in SH coordinates, we have 
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As expected, for v = the above relation is consistent with the radial geodesies of Schwarzschild in standard harmonic 
coordinates [16J ■ 



B. Modified Harmonic Coordinate System 



The SH coordinates are useful for analytical algebraic checks but also contain some gauge-dependent logarithms 
which are less suitable for numerical computations. It has been shown in j25| that such dependences can be transformed 
away by using suitable gauge transformations. The expression for the shift "(^sh^mh)-^" f° r general orbit case has 
been given by Eq. (4.12) of [25| and we have used Eq. (|3.1|1 to obtain the corresponding expression for head-on 
situation. We can write the center-of-mass energy in MH coordinates using the relation 

Emu = E sli + ^sh^mh)^ , (3.10) 

where Egn represents the energy E in SH coordinates and is given by Eq. (|3.6p . "<5(mh-^sh)-E"' for head-on case reads 

~Gm 



J (SH^MH) 



E = 
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log ( ^ I 



■0(8). 



(3.H) 



Case I: Infall from a finite distance 



It is evident from the above that using Eqs. (|3.6[) and (|3.11j) in Eq. f|3 . 10[) we can write the expression for conserved 
energy E in MH coordinates. At the initial separation energy in MH coordinates reads 



Emh(z = Zi) 
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(3.12) 



By equating Eq. fl3.12p and the expression for center-of-mass energy in MH coordinates and then inverting the resultant 
expression, one can obtain the expression for z in MH coordinates. For brevity in presentation, in what follows, we 
will list only the differences in various expressions in a particular coordinate system from their SH values. By adding 
the difference to the SH expression the corresponding expression in the relevant coordinate can be computed. In 
particular in this case, for z we have, 



where, 
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(3.14) 



Though to avoid heavy notation we write z and z, beware that in this subsection they correspond to zmh and zmh 
respectively and in the next subsection to zadm and zadm respectively. 



2. Case II: Infall from infinity 



Once again, by inserting s = z/zi back in Eq. (|3.14p and taking the limit when z, -> oo we obtain the expression 
for '^(sh^mh)^" as 
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z = -V2c 7 7 / 2 
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C. ADM coordinate System 

Finally, in this section we provide the expressions for the conserved energy in ADM coordinates. Like MH coordinate 
systems the ADM coordinate system is also free from logarithms appearing in 3PN expressions of EOM or source 
multipolc moments when standard harmonic coordinate system is used. We can write the center-of-mass energy in 
ADM coordinates using the relation 



-Eadm — Esh + 8 



(SH— >ADM) 



E. 



(3.16) 



where Egn is given by Eq. (|3.6j) and for "<5(sh->adm)^" m head-on situation (see Appendix lAl for its computation) 
we have 
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Case I: Infall from a finite distance 



Eadm(z = z,i) =-fic 7i { 1 - -7, + 
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(3.18) 



For z in ADM coordinate we have, 



ZAVM = ZSR + 0( S H^ADM)Z i 



(3.19) 



where, 
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(3.20) 



2. Case II: Infall from infinity 



The expression for "<5(sh^adm)z" can be written by inserting s = zj ' Zi back in Eq. (|3 . 20[) and taking the limit when 

Zi — > 00 as 



'(SH— >ADM) z 



V2 



1 v 

~ 4 

101959 63 



<V7S I 3-4)7 



5 

16 



3360 64 



16 



--TT^ + lllog - 



(3.21) 



10 



D. Inputs for the computation of the hereditary part 

It is evident from Eqs. (|2 . 14|) - (|2 . 1 5[) that all integrals need to be evaluated with just Newtonian order accuracy 
except the one in the first term of Eq. (|2.14[) which needs to be computed with 1PN accuracy and hence in this 
section we provide 1PN accurate inputs which will be required for the computation of hereditary part of the energy 
flux. Since at 1PN order the expressions for all desired inputs [e.g. source moments, trajectory of the problem and 
relation connecting ADM mass to total mass (m = mi + m^)} are the same in all the three coordinate systems we 
need not give these inputs in different coordinate systems. 



1. Case I: Infall from a finite distance 
Equation (|3.8[) gives the expression for z which in the 1PN limit can be expressed as 

s/T~. s l 



2Gm 



G\ 



- - + - v + s -v 



1 9 



Solving Eq. p.22j) we get the 1PN trajectory as 



3/2 



i \ 1 Gm u i \ 1 Gm u i \ 



(3.22) 



(3.23) 



where g(s), ho(s) and hi(s) are the following simple linear combinations of the elementary functions fi(s) — \/s\/l — s 
and /2(s) = arcsiny^s, 



g(s) = /i(«)-/ a («), 
h (s) = /i(«) + 9/ a («), 

hi(s) = -i(9/i( S ) + / 2 (s)) 



(3.24a) 
(3.24b) 

(3.24c) 



2. Case II: Infall from infinity 



In the 1PN limit Eq. 



gives 



2Gm 



5 Gm / v 
~ 2 ~c Y z V 2 



By integrating Eq. p. 251) we get the 1PN trajectory of the problem which reads 
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IV. THE MULTIPOLE MOMENTS OF COMPACT BINARY SYSTEMS 



In this section we shall provide the expressions for source multipole moments with an accuracy sufficient for the 
computation of the 3PN accurate energy flux in standard harmonic coordinates. General expressions for these moments 
have been given in [25| for inspiralling compact objects moving in generic orbits in standard harmonic coordinates. 
Since the head-on collision problem has only one direction of motion one can write the expressions for source moments 
for head-on case by imposing the restrictions, Eq. (|3.1[) on the corresponding expression for general orbits in terms of 
z and i. Further, as discussed in (29|, the 2.5PN term for general orbits in SH coordinates (see Eq. (3.1) of (25|) is 
modified by 

16 G 2 m 3 v 2 

SJij = - — [-Pfn^rij) + (3 + 3a- 2/3)n^v^] , (4.1) 
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which for the head-on case reduces to —^§- Gr % v z (1 
head-on case reads as 
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(4.2) 



Note that the quantity zq appearing in above expression is, in our present head-on notation, the constant length 
scale ro appearing in Eq. (|2.6|) and in the relations connecting radiative multipolc moments and source multipolc 
moments. The presence of the other constant z' through some logarithms log(z/z' ) at 3PN order is due to the use 
of standard harmonic coordinates and corresponds to r' in the earlier papers. Later in this paper we shall show that 
alternatively one can use other coordinate systems such as the MH or ADM coordinate system which do not involve 
such logarithms. 

The 2PN mass octupole 1^ for head-on case is given by 



hjk = ~mz 3 Vl-4: 
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The 1PN mass moment, Iiju reads as 



'ijkl 
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The moment, Iijkim which will be needed with Newtonian accuracy is 

hjklm = —m z b Vl - 4 v (1 - 2 v) v n^njn k nin m ) + 0(2) . 



(4.4) 



(4.5) 



12 



Finally wc give the monopole moment W, which is 



W = -vmzz + 0(2) 



(4.6) 



V. INSTANTANEOUS CONTRIBUTIONS IN THE ENERGY FLUX FOR HEAD-ON INFALL 
A. The 3PN instantaneous energy flux in Standard Harmonic Coordinates 



Having source multipole moments given by Eqs. (|4.2[) - (|4.6[) and equations of motion given by Eq. (|3.5[) with the 
desired PN accuracies one can compute the required time derivatives of source moments to get instantaneous contri- 
bution to the far-zone GW energy flux using Eq. (|2.8p . Since the instantaneous contribution to 3PN far-zone energy 
flux for compact binaries moving in general orbits has already been listed in [2o| in terms of SH variables r, f and v, 
we can also directly write down the corresponding expression for instantaneous energy flux in terms of the variables 
z and z for the head-on situation using Eq. (|3.1[) . The form of the 2.5PN terms in the flux for a general 2.5PN gauge 
is discussed in [29| [see Eq.(3.14a) there], and we adapt it for the head-on case. We write the result as, 
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and find that the various PN pieces for the head-on case are given by 
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(5.2a) 
(5.2b) 

(5.2c) 
(5. 2d) 
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(5.2e) 



The dependence of the result (|5.1[) - (|5.2[) on z' is due to our use of the SH coordinate system. We will transform 
away this dependence by making use of a different coordinate system such as MH coordinate system. The presence of 
constant zo is not surprising as it was present in the expression of the mass quadrupole moment and hence appears in 
the final expression for the instantaneous part of the 3PN energy flux. This dependence of the instantaneous terms 
on the constant zq should exactly cancel a similar contribution coming from the tail terms. We explicitly show this 
cancellation in the next section. 

The general expression for the energy flux above for the head-on case takes a simpler form for radial infall from rest. 
In this case the velocity z can be expressed solely in terms of the coordinate z and the initial coordinate separation 
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where it is at rest as shown in Scc. lIIII However, since we are working up to 3PN there is one last element to be taken 
into account before we can proceed. This relates to the infall velocity due to leading gravitational radiation reaction 
that induces a zrr. at 2.5PN i.e. i^.sPN- To evaluate this we can adapt the treatment in [HI to the head-on case. It 
requires only the leading term in E and the GW energy flux T and the infall due to radiation reaction for the finite 
separation case is given by 

dE/dt T 16 G 3 m 3 

ZRR = dEjd~ Z = ~dEjd~z = -l5^" (1 ~ S) ■ (5 ' 3) 

Adding zrr to the z given by Eq. (|3.8|) yields the complete 3PN accurate i that will be employed in the next subsection 
to rewrite the energy flux solely in terms of the variable z. It should be obvious that the form of zr.r is the same in 
all the three coordinate systems that we use in this paper. 

We now have all basic ingredients for the computation of the 3PN GW energy flux from compact objects with 
arbitrary mass ratios falling radially towards each other and can proceed to compute the instantaneous part of the 
energy flux for the head-on infall case. 



Case I: Infall from, a finite distance 



Starting from the 3PN instantaneous contribution to energy flux in SH coordinates [Eq. (|5.ip - (|5.2p ] in terms of the 
variables z and z and substituting the expression for i given by Eq. p.8[) supplemented by zrr given by Eq. (|5.3[) we 
get the final expression for energy flux in terms of just one variable z. The final expression for the instantaneous part 
of energy flux in standard harmonic coordinate reads as 
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2. Case II: Infall from infinity 

As discussed in the previous section, we can see case II as a special case of case I and the expression d£/dt can be 
obtained by inserting s = zj z% in Eq. (|5.4[) and taking the limit when 2, — > 00 . The instantaneous part of energy flux 
in SH coordinates reads, 
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The standard harmonic gauge at 2.5PN corresponds to a = — 1 and fj = 0. 

B. The 3PN instantaneous Energy Flux in Modified Harmonic Coordinates 



As we have pointed out in Sec. MI Bl that one needs to use an alternative coordinate system such as MH coordinate 
system, which is more suitable for numerical computations. In this section we provide the 3PN energy flux expressions 
in MH coordinates. We can write the energy flux T in the MH coordinates by using the relation 



■^MH = ^SH + °~ 



(SH-s-MH) 
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where J-sh is the energy flux in SH coordinates for head-on situation given by Eqs. (|5.1[) - (|5.2I) . The general expression 
for the shift "^(sh^mh)^ 7 " is given by Eq. (6.8) of [25[, which in the head-on situation reduces to the following form 
when we use the restrictions given by Eq. (|3.1|) . 
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Once we have expressions for the energy flux T in MH coordinates, we can compute the final expression for instan- 
taneous part of energy flux in MH coordinates following the procedure adopted in Sec. IV Al 



I. Case I: Inj "all from a finite distance 



Substituting for z in MH coordinates, in the expression for energy flux in MH coordinates, one can obtain the 
expression for the instantaneous part of far-zone radiative energy flux in MH coordinates as a function of the separation 
between the two objects at some instant. Rather than writing the full expression for energy flux in MH coordinates 
we list here the difference to be added to the expression in SH coordinates [Eq. (|5.4[) ] to obtain the corresponding 
expression in MH coordinates. We have, 
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2. Case II: Infall from infinity 

By inserting s = z/zt back in the Eq. (|5.8[) and taking the limit when — > oo, the expression for "<5(sh->mh)^/^'' 
takes the form 
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C. The 3PN instantaneous Energy Flux in ADM Coordinates 



In this section we provide the expressions for instantaneous part of the energy flux in ADM coordinates which 
could be useful for the comparison with the numerical relativity results. One can write the energy flux JF in the ADM 
coordinates as 

•^ADM = -^SH + ^(SH^ADM)-^ j (5.10) 

where Jsh is given by Eqs. (|5.1[) and (|5.2[) and the shift "^(sh-j-adm)^ 7 " f° r the head-on situation is 
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We have made use of Eq. (|3.1|) to get Eq. (|5.11|) from the general expression for "^(sh-^Adm)^ 7 " j given by Eq. (6.11) 

of m. 



Case J: Infall from a finite distance 



Substituting for z in ADM coordinates in the expression for energy flux in ADM coordinates, we obtain the 
expression for the difference u S adm) d£ /dt," which should be added to Eq. (|5.4[) to obtain the instantaneous part 
of the energy flux in ADM coordinates. It reads 
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2. Case II: Infall from infinity 



Inserting s = z/zi back in Eq. (|5.12[) and taking the limit when z% 
u S(sH^ABM)(d£/dty as 



00 we obtain the expression for 
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VI. HEREDITARY CONTRIBUTIONS IN THE FLUX FOR HEAD-ON COLLISION 



In this section we shall compute the hereditary contributions to the GW energy flux at 3PN order given by the 
Eqs. (|2.9p . Eqs. (|2.14p . and (|2.15|) . The first hereditary contribution to the energy flux occurs at 1.5PN order and 
is due to GW tails caused by the interaction between mass quadrupole moment and the ADM mass of the source 
causing the spacetime curvature. This contribution is given by the first term in Eq. (|2.14[) where as the second term 
represents the sub-dominant tail at 2.5PN order caused due to interaction of a higher order multipolc moment with 
the ADM mass of the source. Two cubic order tail terms, given by Eq. (|2.15[) . known as tails-of-tails and tail-squared 
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occur at 3PN order and are caused due to the interaction of tails with ADM mass of the source and interaction of 
tails among themselves. 

It should be evident from Eqs. (|2.14j) and (|2.15[) that the computation of all terms would require only Newtonian 
order inputs except the mass-type quadrupolar tail term- first term in Eq. (|2.14p -which would include 1PN correc- 
tions. Note that the second term appearing in Eq. (|2.14|) and needed to be evaluated with Newtonian accuracy, does 
not contribute to the energy flux for the case of radial infall from infinity. The reason is that this term involves 4th and 
6th derivatives of octupole moment [see Eq. (|4.3I) for the expression] but one can check that, at the Newtonian order 
third and higher order derivatives of octupole moment vanish (after we substitute for the Newtonian z appearing in 
corresponding expressions), i.e. I^k = for n > 2, and hence the second integral would not contribute. In this case 
we only have the first term left in Eq. (|2.14[) which gives a hereditary contribution to GW energy flux at 2.5PN order. 
Computation of this term will require the knowledge of 1PN accurate expression for the quadrupolc mass moment 
and 1PN accurate trajectory of the system. In addition to this the 1PN accurate expression for ADM mass would 
also be needed for the computations of tails at 2.5PN order. We have provided the 1PN trajectory in Sec. IIIIDl which 
will be used in computing the hereditary contributions. The second term of Eq. (|2.14p survives for the case of infall 
from a finite initial separation and thus must be taken into account while computing the tails at 2.5PN order for the 
finite initial separation case. As for the instantaneous part we will compute the hereditary contributions as well for 
two different situations, case I- infall from a finite distance and case Il-infall from infinity. 

It is important to note that, at 3PN order unlike the instantaneous part of the flux, the hereditary part is the same 
in all the three coordinate systems - SH, MH and ADM since it involves the inputs which are at most required at 
1PN order and are same in all three coordinate systems. 

In addition to the inputs listed in Sec IIII Dl we also need 1PN accurate expressions for mass quadrupole moment 
and for ADM mass. The mass quadrupole moment in terms of the variables z and z is given by Eq. (|4.2[) . In the 1PN 
limit it reads, 
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The relation between the ADM mass M and total mass m = mi + ?7i2 is given by, 
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where i is given by Eq. 



and is needed to be just Newtonian accurate. 



A. Case I: Infall from a finite distance 



In this case the 1PN accurate expression for quadrupole moment and Newtonian accurate expression for octupole 
moment take the form, 
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lijk = —vmz \]\- Aiyn^rijn k ) . 
The relation between ADM mass and total mass at 1PN order is given as, 
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From the above expression the ADM mass M is independent of z which is consistent with the constancy of M and 
the recognition of the expression of the initial Newtonian energy —Gmvjzi. 
Hereditary terms will involve two integrals, which are given as 



^Quadrupole 
'tail 



£/^M l06 (£ ( "- T) ) + M 



3.5a) 



Octupole 



-'tail 



/ drl\%{r) log (£( 

Ju{Zi) L KZz 



97 
60 



17 
(6.5b) 



Ii 



tail(tail) 



drl[f{r) 



log 2 (£ («-r))+ ^ log (£(«-r)) + 



124627 



44100 



After evaluating the integrals we get 
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where Intl(s), Int20(s), Int21(s), Int30(s), Int31(s), Int40(s). and Int41(s) appearing in Eq. (|6T6al) and Int7(s) 
appearing in Eq. (|6.6b|) are given by 
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and Int5(s) and Int6(s) appearing in Eq. (|6.6cp are given by 
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With Eq. (|6.6[) we can write the various pieces of the hereditary contributions to GW energy flux given by Eq. (|2.9j) 

as 
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Combining Eqs. (|2.9|) . (|6.9|) . (|6.10[) and (|6.11l) now we can write the hereditary contribution at 3PN order involving 
contribution from tails, tail-of-tail and tail squared terms as 
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As we can see the above equation still has some dependence on the arbitrary scale zq at 3PN order. Recall, 
the presence of a logarithmic dependence on zq in the instantaneous contribution to the energy flux at 3PN 
order. The term appearing in the hereditary contribution exactly cancels with similar terms present in instantaneous 
flux expression for energy flux and thus the total flux becomes independent of the arbitrary length scale zq as expected. 



B. Case II: Infall from infinity 



The 1PN accurate expression for quadrupole moment reads 
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It is evident from Eq. (|3.25p . the relation connecting ADM mass M and total mass m [Eq. (|6.2[) ] reduces to 

M = m. (6.14) 

This is consistent with the earlier comment and corresponds to energy vanishing initially. In order to compute the 
hereditary contribution first we need to evaluate the two integrals appearing in Eqs. (|2.14[) and (|2.15p . The integral 
associated with the first term of Eq. (|2.14p is 
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Having all the relevant inputs at the required PN order the value of these integrals read, 
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where the quantity -i/A 1 ) (ii) appearing in Eq. (|6.17b[) is a PolyGamma function whose numerical value is 0.313 25. 
(Of course, formally they correspond to s — > case of the previous section.) 



Using Eq. (|6.17l) in Eqs. (|2.14j) - (|2.15[) we write various pieces of hereditary contribution given by Eq. 
as 
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As mentioned earlier the 2PN accurate energy flux has been given in (l6| which involves the hereditary contribution 
to the energy flux at 1.5PN order. On comparing our results [1.5PN term in Eq. (|6.18[) above with coefficient -5] and 
[Eq. (2.31) of [ill with coefficient -15] for the contribution due to dominant tail we find a mismatch. This apparent 
discre pan cy is a gauge-artifact arising from the difference in our choice of u — t — r/c in contrast to the choice in 
SPW [HI u SPW = t - r/c - (2 G M/c 3 ) log (c 2 r/Gm). We have checked that once we adopt the SPW definition of 
retarded time in harmonic coordinates itspwi our result also leads to the coefficient —15 as in fl(| . This difference 
serves to remind us that the representation of the energy flux in terms of 7 is not gaugeinvariant. 
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Now we can write the total hereditary contribution up to 3PN order to the energy flux as 
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The presence of the arbitrary scale zq in the above expression is similar to the one already noted in Eq. ([6 
will disappear from the final expression for energy flux. 
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12]) and 



VII. THE COMPLETE 3PN ENERGY FLUX FOR HEAD-ON SITUATION 



A. Case I: Infall from finite a distance 



Having computed both the instantaneous and the hereditary contributions to the energy flux at 3PN order for head- 
on situation we are now ready to write the complete 3PN far-zone energy flux due to head-on infall of two compact 
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objects with arbitrary mass ratios. Since the ADM coordinates are independent of gauge-dependent logarithms they 
may be better suited for comparison with numerical relativity results, and we exhibit the complete 3PN accurate 
energy flux expression in these coordinates obtained by adding the hereditary part [Eq. ()6.12j) ] and instantaneous part 
[Eq. (|5.12[) ] of the energy flux. The final result is 
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We can see the final expression for the energy flux [Eq. (|7.ip ] is independent of the arbitrary length scale z . Similarly 
by using Eqs. (|6 . 1 2[) and (|5.4p [ (|5.8p ], one can find the complete 3PN expression for energy flux in SH [MH] coordinates. 
Given the total energy flux as a function of the separation between the two objects at any instant the total energy 
radiated during the infall can be computed as, 
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where Zi and Zf are the initial and final separation between the two objects under head-on infall. Inserting s = z/zi 
and 7 = Gm/c 2 z back in Eq. (|7.1[) and then using it along with i in ADM coordinates in Eq. (|7.2[) one can compute 
the total energy radiated during the radial infall of the two objects from a initial separation Zi to a final separation zj. 
Since Eq. (|7.ip involves some integrals which can only be evaluated numerically, we use the NIntegrate option inbuilt 
in Mathematica to compute the total radiated energy during the process of infall. On the other hand for the case of 
infall from infinity, since we have computed the energy flux as a function of the separation between the two objects 
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in closed form we shall provide 3PN expression for the total energy radiated during the radial infall from z% = oo to 
the hnal separation zt in Sec. lVIFBl however we wish to plot the curves corresponding to the limit Zj = oo with those 
corresponding to the case of infall from a finite distance for comparing the results. 



B. Case II: Infall from infinity 



For this case the complete 3PN expression for energy flux in ADM coordinates can be obtained by adding hereditary 
part (Eq. (|6.2ip ) and instantaneous part (Eq. (|5 . 1 3p ) of energy flux and we have, 
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We can now see the final expression for the energy flux [Eq. (|7.3[) ] is independent of the arbitrary length scale Zq. 
Similarly by employing Eq. (|6.21[) with Eq. (|5.5|) and Eq. (|5.9|) one finds the complete 3PN expression for energy flux 
in SH and MH coordinates, respectively. 

Given total energy flux as a function of the separation between the two objects at any instant the total energy 
radiated during the infall can be computed as 



(7.4) 
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Using the expression for energy flux in ADM coordinates given by Eq. (|7.3j) and z in ADM coordinates given by 
Eq. (|3.2ip in the above we get the 3PN expression for total energy radiated due to head-on infall of two compact 
objects from infinity to a final separation of Zf as 
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where jf = Gm/c 2 Zf. 



VIII. DISCUSSIONS AND CONCLUSION 



Having listed the complete 3PN expressions for the GW energy flux [Eq. (|7.1[) and (|7.3p ] in ADM coordinates, in this 
final section we examine its general behavior as a function of the separation between the two objects under the radial 
infall. Figured] shows the variation of the energy flux, in units of v 1 scaled by a factor c 5 /G = 3.63 x 10 52 joules-sec -1 , 
as a function of the parameter 7 in ADM coordinates (recall 7 = Gm/c 2 z where z is the instantaneous separation 
between the two compact objects falling radially towards each other). Each panel in Fig.[T]shows a comparison between 
the energy flux emitted as a function of the parameter 7 for different initial separations including the limiting case of 
infinite initial separation as well. In each panel curves corresponding to different initial separations (characterized by 
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FIG. 1: Energy flux in ADM coordinates, in units of v 2 (where v is the symmetric mass ratio of the binary), as a function 
of the parameter 7 = Gm/c? z for the head-on situation, for four different initial separations characterized by the parameter 
7; = Gm/c 2 7, = 0.05, 7, = 0.02, 7^ = 0.01 and 7* = 0.0 (infinite initial separation limit) which correspond to the 
situations when the initial separation z; between the two objects is 20Gm/c 2 , 50Gm/c 2 , 100Gm/c 2 and 00 respectively. 
Curves in the top panels correspond to the value of v = (test-body limit) while those in the bottom panels correspond to 
the value of v = 0.25 (equal-mass case). Left, middle and right panels in both top and bottom panels correspond to the 
2PN, 2.5PN, and 3PN accurate expressions for energy flux. The values given on the y-axis have been scaled by the factor 
c 5 /G = 3.63 x 10 52 joules -sec . The labels on the x-axis and alternative x-axis corresponds to inverse of the parameter 7 
(which is the separation between the objects under radial infall at any instant in units of Gm/c 2 ) and the values of the PN 
parameter 7, respectively. 



the parameter 7, = Gm/c 2 Zi) have been plotted for 7, = 0.05, 0.02, 0.01, and 0.0 and correspond to the situation when 
the initial separation Zi between the two objects is 20 Gm/c 2 , 50 Gm/c 2 , 100 Gm/c 2 and 00 respectively. Curves 
in the top panels correspond to v = (test-body limit) while those in the bottom panels correspond to v = 0.25 
(equal- mass case). It is obvious from the figure that the curves in each panel approach each other with increasing 7 
i.e. when the separation between the two objects decreases. This feature can be understood by recalling that since 
s = z/ Zi = 7i/7, for a fixed Zi, the finite-separation corrections in powers of s become progressively less important as 
the bodies approach each other (small z). The finite-separation effects, important when the objects are far apart, are 
less significant at closer separation and the curves for the energy flux approach each other. 

Figure [1] also compares the results that would be obtained using the 2PN, 2.5PN and 3PN accurate expression for 
the energy flux and thus illustrates the improvements arising from a more accurate expression for the energy flux. It 
is clear from Fig. [T] that the energy flux emitted at any instant monotonically increases as the separation between the 
objects under the infall decreases (with increasing 7) as generally expected. However from Fig. [T] we see that after 
a certain maximum value of the parameter 7 in the 2PN and 2.5PN cases the curves show a turnover and start to 
decrease. This is an indication of the fact that the PN approximation is no longer valid beyond this value of 7. It 
should be noted that the value of 7 where this happens depends upon the choice of the initial separation between the 
two objects, the PN accuracy of the expression for the energy flux and the symmetric mass ratio of the binary. 

Finally, Fig. [2] shows the total energy radiated [as discussed in the previous section for the finite initial separation 
case it has to be computed numerically using Eq. (|7.2[) but for the infinite initial separation case it is given by Eq. (17.41) ] 
during a radial infall from initial separation Zi (characterized by the parameter 7$ = Gm/c 2 zi) to a final separation 
Zf (characterized by the parameter 7/ = Gm/c 2 z). Similar to Fig. Q]m Fig. [2] we study the effect of using different 
PN-accurate expressions for energy flux and also the effect of assuming different initial separations in the problem. 
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FIG. 2: Similar to Fig.[T]but for total energy radiated in units of v 2 c 2 m, during the head-on infall of two compact objects from 
a initial separation Zi in ADM coordinates (related to the parameter 74) to a final separation of Zf (corresponding parameter 
7/ = Gm/c 2 z f ). 



It is evident from each panel of the Fig. [5] that as 7/ (zf) increases (decreases) all curves approach each other which 
implies that most of the contribution comes from the late stages of the infall. It is evident from the plots in Fig. [5] 
that only beyond a certain minimum separation between the two objects (under the infall) the estimates of energy 
radiated using PN expressions are valid. The 2PN, 2.5PN and 3PN estimates of the total energy radiated during the 
radial infall (from infinity) of two equal mass compact objects is of the order of 2.2 x 10 -5 , 4.3 x 10 _5 and7.4 x 1CP 5 
respectively. In the test particle limit The corresponding 2PN, 2.5PN and 3PN accurate results for total energy 
radiated in the test particle limit are of the order of 1.4 x 10~ 5 , 3.1 x 10 _5 and 8.5 x 10~ 5 respectively. Unlike the 2PN 
and 2.5PN cases where the breakdown of the PN approximation is explicit in the turnover, the 3PN approximation 
does not show any sharp turnover. As a consequence the value quoted for the maximum energy radiated in the 
3PN case is a bit arbitrary and corresponds to the value at the point where the 2.5PN approximation breaks down. 
From the Fig. [2] one can infer that the energy radiated in the process of head-on infall for the hnitc separation cases 
(7,; = 0.05,0.02,0.01) is of the same order as in infinite initial separation case (7^ = 0). It is evident from the above 
discussion that the 3PN estimates of the peak luminosities and the energy loss in form gravitational radiation during 
the infall between the initial (zt) and a final point (zf) will not only be more than the estimates of the same using 
a less accurate expressions (2PN and 2.5PN accurate) but also they are valid till later stages of the infall and thus 
allows one to compare the results obtained using numerical relativity within the range in which PN approximations 
are valid. 



Appendix A: Calculation of (5(sh— >adm) ^ 

General expression for energy E in CM frame associated with SH coordinate system is given in terms of the natural 
variables; r, v and r [28j]. Noticing this functional dependence and the fact that it is a scalar quantity we expect that 
under a transformation (r' — > r + 5r, v' — > v + Sv, r' — > r + Sr ) this would transform in CM frame as 



E' = E + 5E 



(Al) 



25 



Or equivalently for transformations between SH and ADM coordinate systems, 

EaDM = -EsH + ^(SH^ADM)-S , 

where, 



(A2) 



dE 



dE 



.dE 



^(SH-s-Adm)-E = 8r— + 5v— + 5r — 



(A3) 



The shifts in the variables r, v and f connecting ADM and SH coordinates are given by Eq. (6.10) of |25| an d the 
expression for CM energy Ssh for general orbits is given by Eq. (4.8) of (28j. Having all inputs we now can write the 
shift (5(adm^sh)£' for general orbits which reads as 
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It is easy to see that when restrictions given by Eq. (|3.1[) arc imposed, the above expression reduces to the form given 
by Eq. (I3~T7D . 
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